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Abstract 

We study the gravitational vacuum star (gravastar) configuration as proposed 
by [1] in a model where the interior de Sitter spacetime segment is continuously 
extended to the exterior Schwarzschild spacetime. The multilayered structure 
of [2] - [4] is replaced by a continuous stress-energy tensor at the price of intro- 
ducing anisotropy in the (fluid) model of the gravastar. Either with an ansatz 
for the equation of state connecting the radial pr and tangential pt pressure or 
with a calculated equation of state with non-homogeneous energy/fluid density, 
solutions are obtained which in all aspects satisfy the conditions expected for 
an anisotropic gravastar [IJ. Certain energy conditions have been shown to be 
obeyed and a polytropic equation of state has been derived. Stability of the 
solution with respect to possible axial perturbation is shown to hold. 
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INTRODUCTION 



Recently it has been proposed by several authors that objects other than black holes 
could be formed by gravitational collapse of a massive star. Black hole horizons intro- 
duce number of theoretical problems and a consensus of solutions of those problems 
has not yet been reached [2], [5]. 
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One proposal, which was initiated recently by Mazur and Mottola (M-M) [2]-[l] is 
the so called "gravastar" . In this scenario, quantum vacuum fluctuations are expected 
to play a non-trivial role in the collapse dynamics. A phase transition is believed 
to occur yielding a repulsive de Sitter core which helps balance the collapsing object 
thus preventing horizon (and singularity) formation [S], [S]- is expected that this 
transition occurs very close to the limit 2m{r)/r = 1 so that, to an outside observer, 
it would be very difficult to distinguish the gravastar from a true black hole. Since 
this proposal, different versions of the original gravastar model have appeared with all 
variety of ingredients. For example, it has been shown how a gravastar structure can 
form from a Born-Infeld scalar field such as that predicted by low energy string theory 
[7] or supported from non-linear electrodynamics [8]. In the multi-layer structure of 
the M-M model with de Sitter core and asymptotical Schwarzschild outside region, 
additional features have been added. In [9| gravastar solutions have been studied in a 
generalized (Reissner-Nordstrom) exterior and solutions of the model stemming from 
the original Mazur-Mottola model have been analyzed. Gravastar type solutions in 
the context of solutions to Einstein's equations with tube-like cores have also been 
recently considred [10]. Pioneering in-depth studies of spherically symmetric systems 
with deSitter asymptotics may be found in [11] - [T3]. More recently several papers 
have appeared discussing limits on gravastars and how to distinguish them from black 
holes US], US]. 

It has recently been shown p!] that the gravastar configuration (see the next section 
for the elaborated physical model of it) has to have anisotropic pressures which in ad- 
dition should obey some of the energy conditions of General Relativity. In this context 
one would ideally construct such a spherically symmetric model of an anisotropic fluid 
with a corresponding equation of state. The model, by the definition of the gravastar, 
should not possess a horizon, it should to be stable and its (anisotropic) pressures and 
density ideally would not violate energy conditions. Certainly this last requirement is 
somewhat relaxed in the sense that configurations constructed in this way will violate 
some of the (usual) energy conditions from its very initial definition. Full reviews of 
gravastar models may be found in [T7] and [T5] . 

In this paper we present solutions for the gravastar as proposed by Cattoen, Faber 
and Visser (CFV) [1]. In their paper an attractive sketch (see Fig{T] drawn here for 
the convenience of the reader) is given as a guide to all the future gravastar model 
builders. In the next section a general description of gravastars as anisotropic fluid 
spheres is presented. 

In section 3 we describe the evolution of the anisotropic gravastar model and one of 
the solutions (model 1) which reproduces the CFV sketch of the anisotropic gravastar. 
The value of the surface redshift is calculated while a more elaborate analysis of this 
important quantity is relegated to a future work. Section 4 is devoted to the solution 
(model 2) which starts from the construction of an equation of state and consistently 
solving the remaining generalized Tolman-Oppenheimer-Volkoff equation. Stability of 
this solution with respect to axial perturbations is calculated. Finally, we conclude 
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Figure 1: Sketch of the gravastar as proposed by Cattoen, Faber and Visser [Ij. 



with comments and possible extensions of the anisotropic gravastar solutions. 
2 FROM ISOTROPY TO ANISOTROPY 

Einstein's field equations, being exceedingly complicated because of their nonlinear 
character, have one most important closed form solution outside a spherical star of 
total mass M, namely the Schwarzschild metric. Its line element in curvature coordi- 
nates is 

c/s^ = - 1 \dt^ + — - + dn^ , (1) 

\ r J (1 -2M/r) ' ^ ' 

where dO? ■= dO"^ + sin^ 9 d(f)^. 

For the interior of the star one has to choose a "physically reasonable" stress-energy 
tensor. One attractive possibility is to use the perfect fluid model of matter where 

n = {p + p)u^u,+pS^,, (2) 

with p and p the energy density and pressure respectively in the co-moving frame of 
the fluid, and u'^ being the fluid four-velocity. In the static case, which is studied here, 
the interior metric may be written as 

ds^ = -e<') de + e^(") dr^ + dn\ (3) 
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The Einstein equations, = SttT'^ , give a system of equations: 

(- + \] -\ = 87rp, (4a) 



/y /y'2i j ty'2i 

e-' Q." - IaV + + ^) = Snp , (4b) 

e-" (--\)+\ = 8np. (4c) 

These are supplemented with the conservation law T^.^ = 0, which in this case yields 
only one non-trivial equation: 
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T'\ + - 1 + -Z/' Ti - -A'T° - -T^ = . (5) 
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Elimination of the function vir) from the above under-determined system leads to 
a convenient form of the conservation equation i.e. the Tolman-Oppenheimer-Volkov 
(TOV) [19], [20] equation: 

dp{r) [p(r) +p(r)][m(r) + 47rp(r)r^] 
r^[l — 2m(r)/r] ' 

with 

m(r) := -47r / T\{f)r'^df. (7) 



One may also specify an equation of state which relates pressure and density. If, for 
simplicity, we adopt at this stage a constant density profile function (with built-in 
boundary conditions) 

^^") = {o° for'/Zr 

then it turns out to be an oversimplification which leads to analytic integration of 
([6]) and the corresponding field equations but pressure and density do not obey the 
energy conditions which are analogous to the requirement of mass positivity in the 
Newtonian mechanics. 

With the isotropic fluid and the above homogeneous energy density static solutions 
are allowed for objects with a restricted total mass M to radius R ratio i.e. 2M/R < 
8/9 [21]. 



2.1 ANISOTROPY 



The idea of anisotropy in the spherically symmetric geometry was perhaps first intro- 
duced by G. Lemaitre [22] and suggested by Einstein (as quoted in Ref. [22]). The 
limiting case of p,. ^ is mentioned there and the remaining transversal pressure was 
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said to be enough to support a (stable) sphere (see also [23]). Further development 
has brought different refinements of the original anisotropy notion (see the papers [21] 
- [22] for studies and further references). 

The perfect fluid requires that the pressure in the interior of a star be isotropic, 
leading to calculations of isotropic polytropes for descriptions of objects like white 
dwarfs or neutron stars. Another option giving more freedom to the equation of state 
within the spherical symmetry is the introduction of the stress-energy tensor which 
is anisotropic in its principal pressures. The anisotropy is sometimes (spontaneously) 
produced by extending the notion of a (perfect) fluid to phenomenological models 
including eg. scalar fields, or strongly interacting matter, although it is not known 
how large this anisotropy may be in realistic models. 

The stress-energy tensor for an anisotropic matter/fluid distribution is given by 

n = (P + Pt) U^U, + Vt K + {Pr - Pt) S^S, , (9) 

where Pr is the radial pressure in the co-moving frame and (again) due to the spher- 
ical symmetry, the angular components are identified and are denoted as transversal 
pressures, pf The vector is orthogonal to the fluid four- velocity {s^u^ = 0). 

In the gravastar model, the core interior is assumed to be given by a de Sitter 
solution so the appropriate pressure/density ratio value equal to minus one should 
be implemented as an initial condition for the density profile function, as well as a 
corrector to an equation of state connecting the pressure and density. 

For the moment we assume constant energy density po. In addition to the (energy) 
density in a prescribed form, a relation between the radial pr and the tangential pt 
pressures should be given. An ansatz for the anisotropy measure, 

A := ^ , (10) 

will be used following the hints given in [1]. Bounds on the anisotropy measure are 
calculated and are expressed in terms of the "compactness" 2m{r)/r, so our ansatz 
has the following form: 

a^Mr) 

12 r ^ ^ 

The constant a^/12 will simplify the (numerical) calculation and it is a measure of 
anisotropy for this version of the gravastar model. The TOV equation now assumes 
the following form: 



r'^Pr 



A-KT 
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a Po 



3[P0 +Pr(r)][po + 'iPr{r)] 



1 - 8r%po/3 

with the initial condition PriS^) = po, requiring also that Pr{0) = pt{0). 
The radial pressure Pr{r) is given by 



(12) 



Pr{r) = y 



-1 +a^/l - p(r) -1 + 



2(a-2) 



a -2 + {a + 2) exp[a(-l + ^1 - p(r))] 



13) 
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where the "compactness" is /i(r) = 2m{r)/r = Snr^po/S, while the transversal pressure 
can be obtained through (fTTj) . and is given by 



.1+ 



- 2 + (a + 2) exp[a(-l + ^1 - /i(r))] 



+ /ia' V . (14) 



For all a > 0, both pressures start at — po in the centre of the gravastar. At /i = 1 (if 
one would allow such a density profile) the values are Pr = — po/3, Pt = ^ (a^ ~ !)• 
Only for values a > 4.11172, does the radial pressure reach positive values before 
returning back into the negative pressure region. It is also worth noting that this 
simple model does not lead to an equation of state, since the energy density does not 
change within the star. 

As previously mentioned, this oversimplified model with the constant (energy) den- 
sity function as given above cannot reproduce the sketched Figure [H i.e. it cannot 
provide a gravastar in the proposed model without the surface layer structure as re- 
quired by junction conditions. The most obvious extension is to use perhaps the same 
ansatz for the equation of state connecting the behaviour of radial and tangential pres- 
sures as given in ( fTTj) . For the density one takes some non- homogeneous distribution 
which will give a more complex dependence of Pr{r) and lead to a possible solution 
(see the model 1 below). 

Certainly the ansatz could be replaced by a calculated equation of state and intro- 
duce again an inhomogeneous density distribution (see the model 2 below). On these 
two possibilities the new results of this paper are based. 

2.2 ENERGY CONDITIONS 

In order to have a reasonably realistic radial pressure component we demand the (most 
natural) weak energy condition (WEC) be satisfied everywhere: 

p(r) > along with p(r) +p^(r) > and p(r) + Pt{r) > . (15) 

Also, we require 

p(r = 0) > (16) 

and, since stability of the fluid/gravastar requires that p(r) must not increase out- 
wards, 

dp{r) 



dr 

The above WEC obviously implies: 



< . (17) 



p(r) + Prir) > and p(r) + Pt{r) > , (18) 
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the so called null energy condition (NEC). There is also a dominant energy condition 
(DEC) which requires that 



p(r) > and Prir) E [—p{r), +p{r)] and Ptir) E [—p{r), +p{r)], (19) 

and which implies the other two energy conditions. 

Another commonly studied energy condition is the strong energy condition (SEC) 
which states that, for our static system, 

P + Pr > 0, p + pt>0 and p + Pr + 2pt > 0. (20) 

Since gravastars possess a deSitter core, it is not possible to satisfy this energy condi- 
tion. 

Within the normal range of the gravastar it could be reasonable to propose that 
the speed of sound shall not exceed 1 (speed of light). This requirement is reasonable 
to apply in the region where one could expect that the unusual physics govern most 
of the physical processes (i.e. where the matter possesses the least exotic behaviour). 
From the sketch in Figure [T] this is expected to be in the gravastar atmosphere, so 



2 dp 

V. — 



dp 



atm. 



< 1- (21) 



In this region one should be able to derive a polytropic equation of state which 
should be expressed by parameters which are in accordance to the requirement involv- 
ing the speed of sound (see next Section). 



3 MODEL 1: THE ANSATZ AND AN INHOMOGENEOUS DENSITY 
DISTRIBUTION 

In this section the solution of the gravastar non-layered model (i.e. the model sketched 
in Figured] [Ij) is sought by prescribing an equation of state relating the pressures, 
i.e. an equation which connects the radial Pr{r) and tangential pt{r) pressure. An 
improved ansatz of the form 

A = A-^ (22) 

Po 

is assumed, where A is as given in (fTTl) . The (inhomogeneous) density function profiles 
will be chosen to be either of the exponential form or a polynomial form. 

The gravastar is a static configuration with a prescribed behaviour of the (radial 
and tangential) pressure. As sketched in Figure [1] the radial (and tangential) pres- 
sure originates at r = with the initial value conforming with the de Sitter interior 
definition {pr{r = 0) = pt(r = 0) = —p{r = 0)). Investigations concerning a possible 
analytic form of the (energy) density have produced stringent bounds on the allowed 
functional behaviour. 
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In the first approach to the gravastar model we start from the anisotropic TOV 
equation. Following the requirement imposed on the density profile we will be using 
a simply behaving function of the form 

p(r,ro)=poe-(^/'^«)\ (23) 

where tq has an appropriate dimension and n is chosen to be a positive integer. Mo- 
tivations for the exponential form of p may be found in [13], [3^1, where they chose 
n = 3, which simplifies some calculations. The values of the central density, Pq, and 
of the fall-off constant, tq, may not be set arbitrarily, since one must make sure that 
the compactness does not exceed unity. Alternatively, one may require a certain total 
mass of the configuration, and set the upper limit to the compactness within the star. 
Then the values for the parameters po and tq would follow. In this (and the following) 
model we will restrict our (numerical) procedures to a general total mass M and tq 
by keeping the corresponding quantities as ratios. 

When requirements of anisotropy are being imposed we use the anisotropic fluid 
stress-energy tensor ([9]) T'^ = diag(— p, Pr, Pt, Pt)- As before, we have from the equa- 
tions fHal) -(B 

m{r) = Atx I p{r)f'^df (24) 
Jo 

and the generalized TOV equation 

= -[p(r) +Prir)]^^j^^^^^^ + 2 ^^^^ 

= - P + Vr{r)\ I / + 2p 

r"'[l — 2m[rj/r\ r 

where we have introduced the anisotropy parameter A(r) from above. 

In Figure [2] the (numerically obtained) solutions of the transversal pt and radial Pr 
pressures are depicted when the modified ansatz given in (1221) is used as an input for 
the equation of state. Also the density profile p(r) given in fl23l) (decreasing dashed 
line) as well as po • p(r) = 2po ■ m{r)/r for po = 1 (initially increasing dashed line) is 
presented. 

The radial pressure Pr{r) as well as energy density satisfy all the energy conditions 
assumed to be valid for gravastar solutions. 

Instead of insisting on numerical details we want here to show that the proposed 
gravastar structure could be found also with other energy density profile functions 
which also follow general requirements for a gravastar. By choosing another indepen- 
dent profile we will show that solutions found in this paper are not just artefacts of a 
very specifically chosen energy density profile but a more general feature of prescribed 
initial conditions within the framework of general relativity. 
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Figure 2: The gravastar model with the energy density profile p = po exp[— (r/rg)"] and 
anisotropy {pt —Pr)/p = (p/Po) /^/12: radial (lower solid line) and transversal (upper 
solid line) pressures, energy density (dashed line) and the compactness (dotted line). M is 
the total mass of the configuration. In this example the parameters are n = 3, maximal 
compactness within the gravastar pmax. = 0.822, and a = 6.69. M = 1 here so, for a body 
of mass 10 km, distances are measured in units of 10 km. 

3.1 MODEL lb 

As the next example of the gravastar model construction we will use the density profile 
function of the form 



and with the ansatz given in (iTTl) . Again the values of the central density, po, and of R 
are adjusted to make sure that the compactness does not exceed unity. As an example 
we used n = 4 in fl26|) to construct a gravastar of total mass M = 1 with maximal 
compactness within the gravstar /^max = 0.93. For these parameters we obtained 
Po = 0.037 and R = 2.31, and the compactness at the surface /i(-R) = 0.864. If M is 
measured in decametres, the central density in this model is approximately 4 x 10^^ 
kg/Dm^, approximately 1000 times neutron star density. The gravitational repulsion 
of the near deSitter type core region is supporting such densities. 

By adjusting the value of a we could obtain a solution in which both pressures 
vanish at r = R, thus obtaining a gravastar. The appropriate value was a = 6.10, 
while for lower/upper values the pressures end up at negative/positive values at the 
surface of the gravastar, such solutions were not considered any further (although the 
tangential pressure and energy density need not vanish at the surface). The pressures, 
the density profile and the compactness of this model are shown in Fig. [31 




n>2 , 



(26) 
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Figure 3: The gravastar model with the energy density profile p = pq{1 — {r/R)'^) and 
anisotropy {pt —Pr)/p = oP' (p/po) ^/12: radial (lower solid line) and transversal (upper 
solid line) pressures, energy density (dashed line) and the compactness (dotted line). M is 
the total mass of the configuration. In this example the parameters are n = 4, maximal 
compactness within the gravastar ^max = 0.889, and a = 6.10. Surface compactness is 
^(R) = 0.864. 



3.2 THE EQUATION OF STATE AND SURFACE REDSHIFT 

Apart from the energy condition relations among the components of T'^, an equation 
of state for the matter described by should exist. An equation of state could 
be either given in the sense that it is prescribed as an intial condition under which 
the Einstein field equations should be solved or derived as a posteriori calculated law 
which should be in accordance with the expected characteristics of a physical system 
described by the T'^ content. 

The dependence of the (radial) pressure on the energy density is usually, in the 
perfect fluid regime, expressed in a mathematical form 

p = Kp^+y^p = Kp' (27) 

which represents the polytrope, and Up is a polytropic index [37]) . 

According to the graph of the radial pressure vs. density (in the mirrored presenta- 
tion - see FigureHj) it is clear that the equation of polytrope could/should be calculated 
only in the "atmosphere" of the gravastar (see Figured]), where radial pressure and 
density follow the usual behavior of a more standard fluid and the equation of state 
is assumed to relate energy density to radial pressure, as is commonly employed when 
dealing with anisotropic fluids. In the polytropic regime the speed of sound respects 
the bound dp/ dp < 1. The specific models constructed here are designed to illustrate 
the maximum compactness allowed by the speed of sound causality condition fl2T]) . 
The parameters can easily be adjusted to produce models where the speed of sound 
is much less than unity. 
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One method of determining a star's properties is the measurement of its surface 
redshift, Z, of spectral hues produced in the star's photosphere. It is defined by the 
fractional change (AA) between observed (Aq) and emitted (Ae) wavelength compared 
to emitted wavelength: 

and, according to our notation, its surface value becomes 

Zj^ = e"^(«)/2 - 1 . (29) 

In the context of the presented model calculations it is important to recall that 
anisotropy affects the (surface) redshift so that for the static perfect fluid sphere the 
surface redshift is not larger then Zs = 2 [21] , [28] whereas in [28] , [31] , and references 
therein, the maximal surface redshift for anisotropic sphere is found to be 3.842. The 
model la (see Fig. 2) provides the value Zg = 1.230 (with maximum compactness 
A^max = 0.822). The model lb (Fig. 3), with /imax = 0.889 and surface compactness 
/is = fi{R) = 0.864 gives the value Zg = 1.712. Both values are well within expected 
values for a realistic (neutron star) object. 




Figure 4: The equation of state resulting from the solution shown in Fig. [2 radial pressure 
(thick solid line) and its derivative with respect to the energy density (thick dashed line). 
The polytropic fit computed in the range indicated by the marker and extrapolated over 
the whole range of densities: radial pressure (thin solid line) and its derivative (thin dashed 
line). The parameters of the polytropic fit are k = 9.16 and 7 = 1.70. The atmosphere 
corresponds to the left region of the graph. 
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Figure 5: The radial pressure of Fig. [2] (solid line) and the polytropic fit radial pressure 
(dashed line) vs. r. 




Figure 6: The equation of state resulting from the solution shown in Fig. [3) radial pressure 
(thick solid line) and its derivative with respect to the energy density (thick dashed line). 
The polytropic fit computed in the range indicated by the marker and extrapolated over 
the whole range of densities: radial pressure (thin solid line) and its derivative (thin dashed 
line). The parameters of the polytropic fit are k = 38.0 and 7 = 2.18. The atmosphere 
corresponds to the left region of the graph. 
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Figure 7: The radial pressure of Fig. [3] (solid line) and the polytropic fit radial pressure 
(dashed line) vs. r. 
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4 MODEL 2 



Here we select a density profile of the form (1251) . Specifically, for numerical calcula- 
tions, we chose 

Po = Y^2^ ro = (12)V3M, M = 6000m. 

The energy density is related to the radial pressure via a Mbonye-Kazanas (MK) 
equation of state [36] : 



s-is + l) 



l/n 

I I 1 1 I I \ I I II I I \ 

Pr{p) 



The MK equation of state possesses several desirable features (with appropriately 
chosen parameters), namely: the speed of sound is less than one in the atmosphere 
and the WEC and DEC are satisfied. The model considered here utilizes (|30l) with 
the parameters m = 2 and n = 1. The tangential pressure will then be given by the 
anisotropic TOV equation ( l25l) . 

In Figure [8] we plot Pr/p (solid) and pt/p (dashed) in order to study the DEC. 
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Figure 8: The DEC with radial pressure Pr{r)/ p{r) and with tangential pressure Pt{r) / p{r) . 
The parameters are: s = 2.2135 and total mass M = 6 km. 



From Figure [H] we see that there exists DEC violation with respect to the tangential 
pressure. This violation occurs in the crust of the gravastar, where the physics is 
expected to be "exotic" and generally cannot be avoided as noted in [I]. Note that in 
a model constructed via this method, the DEC violation is minimal and the DEC can 
be respected in the outer layer of the star. 

The expression for the anisotropy measure A (which is now calculated and not 
prescribed as in model 1), although in closed form, is not simple. Instead, for the 
specific model considered here, we plot it in Figure [H Notice from this figure that 
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A vanishes in several locations. It should be possible at these points to patch the 
solution to a perfect fluid yielding a regular star-like structure in the outer region. 

Following the calculation in the previous model, the specific model constructed 
here possesses a surface redhift of = 0.4 . 



Figure 9: The anisotropy parameter, A = {pt{r) — Prir))/ p{r), for model 2. 



4.1 AXIAL STABILITY 

Here we analyze stability of this model against axial perturbations. Stability of gravas- 
tar models with thin-shells against spherical perturbations has been analyzed in several 
papers [SB] , [3S] . The model of Lobo [3S] considers the interesting scenario where the 
gravastar like object is formed from the gravitational condensation of dark energy, 
believed to be responsible for the current acceleration of the universe. The physical 
motivation for considering dark energy stars may be found in ^0]. 

We concentrate here on the issue of axial perturbations, which allow for pertur- 
bative rotations. This is useful as, in a realistic astrophysical collapse scenario, one 
would expect the star to possess some amount of angular momentum. For details 
the interested reader is referred to the book by Chandrasekhar [H] and the paper by 
Dymnikova and Galaktionov [53] . Before proceeding we shall establish some notation. 
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(32) 
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(these will be used below). 

A more general line element may be written as 



-e" dt^ + dr^ + dO'^ + sin^ 9[d(t)-ujdt- dr - dOf 



(33) 



By comparison with ([3]), the unperturbed metric has u = q2 = = 0. Axial per- 
turbations correspond to these becoming non-zero. Our analysis is to linear order in 
these quantities. 

We study the form of the equations asymptotically, in the region where = —e~^, 
within a time rescaling. With some restrictions, this will be sufficient. A full stability 
analysis, which is general and does not rely on our assumptions may be found in |16] . 

We write 

cj(r, e, t) = cD(r, 5)e*"*, (34) 

and similarly for q2 and q^. 

For time-harmonic perturbations, it can be shown that the system governing the 
perturbations can be reduced to a single second-order differential equation, which 
can be solved by separating the variables r and 6 [33]. In brief, the perturbed field 
equations give a relation between u and q2 and gs, allowing the elimination of u. The 
quantity Q := e''r^sin'^6' {d0q2 — drqs) is written as Q = R{r)Q{9) and the resulting 
radial equation is, asymptotically. 



2 u 

r e 



d 
dr 



dR 

J.2 ^j. 



,e''R 



+ a^R = 0, 



(35) 



where fif is the eigenvalue of the angular equation which can take on the values 
= {l + 2){l - 1) for / = 2, 3, ... . 
We can make a change of coordinates, 

e-"" dr , 

so that ( l35l) reduces to a Schrodinger type equation: 

d^ 



drt 



- Viir) 



(36) 



where 



and 



ZAr] 



Ri{r,) 



Viir) 



{jjf + 2e" - e^ru') 



(37) 



The full potential, valid throughout the domain, may be found in [16j where a full 
stability analysis is performed. An equation of the form (136|) results even away from 
the end points [32] . 
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The allowed values of a are therefore the eigenvalues of (!36l) . If all the eigen- 
values are positive, then a will always be real. Hence the time dependence of the 
perturbations will oscillate but not grow. These perturbations are then stable. 

The eigenvalues, o";, are determined by the potential Vi{r). In fact, knowing the 
asymptotic behaviour of the potential is sufficient for time-harmonic perturbations as 
observed in [13]. The eigenvalues will be positive if Vi{r) is positive and the asymptotic 
behaviour near zero and as r — > oo both go as r~^. 

The positivity of the full V;(r), as given in [16j, is guaranteed if there is no DEC 
violation throughout the star or, at the very least, if the DEC violation is not too 
large. 

Regarding the asymptotics, we consider the limit as r ^ which, for gravastars 
with de Sitter centers gives 

limz/(r) = and limz/'(r) = 0. 

r— >0 1 — >0 

These conditions lead to 

V. - ^ (38) 

as r — i> 0. In the far region, r oo, we need to show that r(r) possesses finite limit. 
For large enough r it can be seen (from Figure [8]) that r — 2m(r) >> 1 and that 
Pri''^) < p{r) hence 

h{r) < 167rrV(r). (39) 

Since the integral of (l39i) has finite limit, r(r) will also have finite limit and therefore 
z/(r) also possesses finite limit. In addition. rh{r) ^ in the far zone so that rz/'(r) 
vanishes. Thus, for very large r, Vi{r) has behaviour. This completes the argument 
for axial stability. 



5 CONCLUDING REMARKS 

In this paper we have constructed several models of gravastars which do not employ 
thin shells. These models obey the conditions posed in [Ij . The "stellar construction" 
is based on two methods. One method employs the prescription of the energy density 
and anisotropy parameter, where the system can be shown to fit a polytropic equation 
of state in the stellar atmosphere. In the second method, a presciption of the energy 
density is imposed, and then the radial pressure is defined via a reasonable equation 
of state. A general model of the latter type has been shown to be stable under axial 
peturbations. It is expected that similar models, utilizing the former method, would 
also exhibit such stability. In closing we should mention that the models presented here 
are proof-of-concept models. That is, they are constructed to show that gravastars, 
under the requirements set out by previous authors (for example [Tj) can be explicitly 
achieved within the framework of general relativity theory. It would be of interest to 
study whether the gravitational collapse of a heavenly body, with reasonable initial 
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conditions and undergoing the phase transition proposed in [2]- [5] at late times, will 
always form an object without an event horizon. 



ACKNOWLEDGEMENTS 

D. H. is thankful to Simon Fraser University for their hospitality and support under 
which a part of this work was carried out. This work was partially supported by the 
Ministry of Science of The Republic of Croatia under the project number ZP0036- 
038 (D. H.) and the Natural Sciences and Engineering Research Council of Canada 
(K. S. v.). We would also like to thank the anonymous referees for their comments. 

References 

C. Cattoen, T. Faber, and M. Visser, Class. Quantum Grav. 22, 4189 (2005) 



P. O. Mazur and E. Mottola, arXiv.org gr-qc/0109035 



P. O. Mazur and E. Mottola, Proceedings of the Sixth Workshop on Quantum Field The- 
ory Under the Influence of External Conditions, Oklahoma, Rinton Press, Princeton, 
N.J. (2003), arXiv.org: gr-qc/0405111 

P. O. Mazur and E. Mottola, Proc. Nat. Acad. Sci. Ill, 9545 (2004) 

G. Chapline, E. Hohlfield , R. B. Laughlin, and D. I. Santiago, Int. J. Mod. Phys. A18, 
3587 (2003) 

E. B. Gliner, Sov. Phys. JETP, 22, 378 (1966) 

N. Bilic, G. B. Tupper, R. D. Viollier, J. Cos. Astrophys. 0602 (2006) 013 

F. S. N. Lobo and A. V. B. Arellano, Class. Quant. Grav. 24, 1069 (2007) 

B. M. N. Carter, B. M. N. Carter, Class. Quant. Grav. 22, 4551 (2005) 

0. B. Zaslavskii, Phys. Lett. B634, 111 (2006) 

1. Dymnikova, Gen. Rel. Grav. 24, 235 (1992) 

I. Dymnikova and B. Soltysek, Gen. Rel. Grav. 30, 1775 (1998) 

I. Dymnikova, Grav. Cosmol. Suppl 8N1, 131 (2001) 

I. Dymnikova, Int. J. Mod. Phys. D12, 1015 (2003) 

A. E. Broderick and R. Narayan, Class. Quant. Grav. 24, 659 (2007) 

C. B. M. H. Chirenti and L. Rezzolla, larXiv:0706.15T3] 
C. Cattoen, arXiv.org: gr-qc/060601l] 



T. Faber, arXiv.org: gr-qc/0607029 



18 



[19] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 374 (1939) 
[20] R. C. Tolman, Phys. Rev. 55, 364 (1939) 
[21] H. A. Buchdahl, Phys. Rev. 116, 1027 (1959) 

[22] G. Lemaitre, Ann. Soc. Sci. Bruxehes, A51, 51 (1933); reprinted in Gen. Rel. Grav. 
29, 641 (1997). 

[23] P. Florides, Proc. Roy. Soc. London A337, 529 (1974). 

[24] R. L. Bowers and E. P. T. Liand, ApJ. 188, 657 (1974) 

[25] G. Rein, Math. Proc. Camb. 115, 559 (1994) 

[26] J. K. Rao, M. Annapurna and M. M. Trivedi, Pramana - Journal of Physics, 54, 215 
(2000) 

[27] L. Herrera and N. O. Santos, Phys. Rept. 286, 53 (1997) 

[28] B. V. Ivanov, Phys.Rev. D65, 104011 (2002) 

[29] T. Harko and M. K. Mak, Annalen Phys. 11, 3 (2002) 

[30] K. Dev and M. Gleiser, Gen. Rel. Grav. 34, 1793 (2002) 

[31] D. E. Barraco, V. H. Hamity and R. J. Gleiser, Phys. Rev. D67, 064003 (2003) 

[32] M. K. Mak and T. Harko, Proc. Roy. Soc. Lond. A459, 393 (2003) 

[33] L. Herrera, A. DiPrisco, J. Martin, J. Ospino, N. O. Santos and O. Troconis, Phys. 
Rev. D69, 084026 (2004) 

[34] H. Hernandez and L. A. Nunez, Can. J. Phys. 82, 29 (2004) 

[35] S. D. Maharaj and M. Chaisi, Math. Methods Appl. Sci. 29, 67 (2006) 

[36] M. Mbonye and D. Kazanas, Phys. Rev. D 72, 024616 (2005) 

[37] S. Weinberg, Gravitation and Cosmology, John Wiley, New York, (1972) 

[38] M. Visser and D. L. Wiltshire, Class. Quant. Grav. 21, 1135 (2004) 

[39] F. S. N. Lobo, Class. Quant. Grav. 23, 1525 (2006) 

[40] G. Chapline, Proceedings of the 22nd Texas Conference on Relativistic Astrophysics, 
California, (2004), arXiv.org: astro-ph/050320 

[41] S. Chandrasekhar, The Mathematical Theory of Black Holes, Oxford University Press, 
Oxford (1992) 

[42] V. Ferrari, Gravitational Waves, Stars and Black Holes, Invited talk, in Black holes 
and Relativistic Stars: A Symposium, University of Chicago, (1997) 

[43] I. Dymnikova and E. Galaktionov, Class. Quant. Grav. 22, 2331 (2005) 



19 



